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I utilize the Caves- Milburn model for continuous position measurements to formulate a broadband 
version of the Standard Quantum Limit (SQL) for monitoring the position of a free mass, and 
illustrate the use of Kalman filtering to recover the SQL for estimating a weak classical force that 
acts on a quantum-mechanical test particle under continuous observation. These derivations are 
intended to clarify the interpretation of SQL's in the context of continuous quantum measurement. 
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Substantial efforts have been devoted to elaborating 
Standard Quantum Limits (SQL's) for both discrete 
and continuous [^J^] measurement of the position of a 
free mass. In large part, the motivation for such investi- 
gations stems from a pressing need to identify any pos- 
sible constraints imposed by the principles of quantum 
measurement on the experimental possibility to detect 
gravitational waves My objectives in this paper 

will be to formulate an SQL appropriate to broadband 
continuous measurements of the position of a free mass 
(expressed in terms of the signal bandwidth and with 
measurement sensitivity expressed in units of length per 
VHz), to consider its proper intepretation, and to demon- 
strate its compatibility with the usual SQL Q for detect- 
ing a weak classical force. The mathematical analysis 
will be adapted to the continuous measurement model of 
Caves and Milburn , which appears to have direct rele- 
vance for concrete experimental scenarios such as atomic 
force microscopy (9| and cavity QED 10 PI]. 

Although the limits I derive will be familiar from previ- 
ous studies of the discrete and narrowband measurement 
scenarios, my emphasis here will be on formulating these 
limits in a manner that is specific to broadband quan- 
tum measurement. In particular, I have found that some 
care needs to be taken in deriving an SQL for detect- 
ing weak forces from the SQL for continuous position 
measurement. The method I use below emulates the 
classical technique of Kalman filtering (ljjll|, and was 
motivated by the general strategy of quantum system 
identification discussed in [Q. Throughout this paper 
I have adopted the time-domain, state-space perspective 
of quantum measurement whose virtues have become ap- 
parent from recent advances in adaptive quantum mea- 
surement Q and quantum feedback jl6| . Such methods 
hold great promise for the challenge of formulating ex- 
perimentally tenable strategies for broadband quantum 
nondemolition measurement. 

The mathematical formalism necessary to treat con- 
tinuous quantum measurements has been developed by 
numerous authors, with the most relevant works for the 
present discussion being ]|,[l3|,|l7]] . While it is not abso- 
lutely essential to go to the continuous limit, doing so 
will allow us to use the convenient notation of stochas- 
tic differential equations (SDE's) ]lS|] . What really mat- 
ters for the discussion at hand is an assumption that the 
timescale associated with measurement interactions and 



readouts is much shorter than any timescale on which we 
wish to undcrtstand the system dynamics. 

In the Caves-Milburn model for continuous quantum- 
mechanical measurements of position |^| there is one sys- 
tem of interest, and an infinite succession of identical 
(and identically-prepared) "meters." Let the meters be 
labelled by an index r. The system is brought into mo- 
mentary contact with the r th meter at time t r — rr, and 
the position operator of each meter is measured sharply 
just after it has interacted with the system. The string of 
measurement results thus generated constitutes a classi- 
cal record of the system evolution. Continuous measure- 
ment is achieved in the limit where the time-interval be- 
tween measurements r goes to zero, with the uncertainty 
of each individual measurement simultaneously going to 
infinity in an appropriate manner. 

Let x,p denote the system position and momentum 
operators, and X ri P r the position and momentum op- 
erators of the r th meter. We adopt the following time- 
dependent interaction Hamiltonian between the system 
and the meters: 
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The delta-function form is chosen to facilitate the limit 
t — > 0, but as discussed above it may be viewed as the 
idealization of any "shrinkable" function with compact 
support. We assume that each meter, just before it in- 
teracts with the system, is prepared in the pure state 
|T r ) with Gaussian wave-function (X r = £ r |£ r )) 



<£ r |T r ) = T(£ r ) = 
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Caves and Milburn have have derived an exact expression 
for the conditional evolution of the system state under 
the following measurement protocol: at each time t r , i. 
couple the system to the r th meter prepared in state (||) , 
ii. evolve the system and meter under the interaction 
Hamiltonian ([!]), Hi. perform a precise measurement of 
X r . If we write p{t r ~) for the system state just before 
the measurement at time t r = rr, the post-measurement 
system state p(t r +) is given by 



p(t r +) =t(f P )p (*,.-) ft (&.), 
T(£ P ) =(C r \e- ixP ^ n |T r >, 
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where £ r corresponds to the actual outcome of the meter 
projection in stage Hi. Recall that x, X r P r are Hermitian 
operators, £ r is a c-number, and | T r ) is a Hilbert-space 
kct. T r is thus a quantum operation density. Qualita- 
tively speaking, the basic effect of (||) is to shift the cen- 
troid x' r = Tr [xp (t r — )] of the system's position-space 
distribution towards the value conveyed by £ r , and to re- 
duce its overall width A' r = Tr [x 2 p(t r — )] — (x' r ) 2 by an 
amount that depends on the ratio A' r /a. Note that this 
is not a projective measurement — the post-measurement 
state p (t r +) depends on both the measurement result £ r 
and the pre-measurement state p (t r —) . 

The operation (Q) maps Gaussian pure states of the 
system to Gaussian pure states, so we can in fact 
parametrize the selective evolution of a free particle (ini- 
tially prepared in a Gaussian pure state) by just four real 

numbers. Chosing x r = (x) ,p r = (p) , A r = 2 ^(A.t) 2 ) , 

and e r = (AxAp — ApAx) /%, Caves and Milburn de- 
rived a set of difference equations for the evolution 
of these quantities with each successive measurement. 
They further showed that Gaussian initial states gen- 
erally converge toward stationary evolutions, in which 
x r and p r evolve stochastically but A r and e r vary 
periodically — during the time intervals between measure- 
ments both widths increase according to the free Hamil- 
tonian, but the effect of each measurement is to reduce 
them by a constant "contraction factor" C such that 
A r = Tr [x 2 p (*r+)] — i x r) is independent of r. The 
form of the evolution equations is especially convenient in 
this stationary regime, so in what follows we shall assume 
that stationarity has been bootstrapped by preparing the 
system in an appropriate initial state. 

In order to derive a broadband SQL for position mea- 
surements, we start with the stationary difference equa- 
tions (3.30a,b) from M: 
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Here unprimed quantities x r , p r describe the system state 



just after the (r — l) th measurement has been made, and 
x' r holds just before the r th measurement. 

We first need to rewrite the difference equations (||) as 
a stochastic differential equation (SDE), by taking r — > 
together with a — > oo such that their product D = <jt 
stays constant. Recall, from equations (3.20,3.21) of ||, 
that (£ r — x' r ) is a Gaussian random variable with zero 
mean and variance aC/2. In the continuous limit it also 
turns out that C — * 1 so we immediately have an 
explicit (Ito) stochastic differential equation: 



where dW is a Wiener increment fll 




dW, 



(6) 



To derive an SDE for the rms position, we begin by 
applying Ito's formula to derive an SDE for x ms = (a; 2 ) , 



d (x 2 ) = — xdt, 
m 



dt ' 



= — (px) 
m 



(7) 



Substituting x,p by integrals of equation (g), and keeping 
in mind that the "functions" dW (t) appearing in the 
expressions for x and p will be identical for any given 
stochastic realization, 



dt ' 



Hence, 
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We now turn to quantify the fundamental measure- 
ment noise inherent in the present scheme. As mentioned 
above, the discrete measurement errors (before taking the 
continuous limit) are a stationary Gaussian process with 
variance Ca/2. Hence, in any given time interval H~ x 



the sample variance for the N 
will be 
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Note that v D thus represents a measurement inaccuracy 
(or "sensitivity," depending on your point of view) in 
units such as meters/ \/Hz. In the sequel we shall regard 
S (t) as the time-domain measurement signal, so that B 
should be interpreted as a low-pass bandwidth. It is im- 
portant to understand that D and B are completely in- 
dependent variables — D parametrizes the system-meter 
coupling strength (a physical quantity), B the degree of 
smoothing applied to the measurement results (a signal- 
processing quantity). 

We can now formulate a condition for backaction- 
induced diffusion to become visible against the funda- 
mental measurement noise: 



t 3 / 2 = ^/DB/2, 
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(11) 



This expression should be taken to predict the overall 
measurement time at which backaction-diffusion 
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Figure 1. (a) Horizontal lines indicate the rms signal 
noise ±y/~(DB /2), curved lines indicate ±x rms (t), and 
the vertical line indicates t* for D — 1.42 x 10~ 20 , B = 
10 7 , and m = 2.22 x 10 -25 (mks units). Two stochastic 
realizations of the measurement signal are also shown, 
generated by direct simulation of equations (19)-(21). (b) 
The nonstationary component of the signal £ r can be 
perfectly subtracted away to yield eta r = £ r — x' r , which 
is a Gaussian noise process with variance aC/2. 

should cause the rms wander x rras of the observed 
particle (with mass m) to equal the point-to-point scat- 
ter S rms in continuous measurements made with mea- 



surement inaccuracy V D and filtered to bandwidth B. 
Viewed as such, I claim that equation ( JTT| ) represents the 
most sensible form for a broadband SQL for continuous 
observation of the position of a free mass. Note that 
the timescale for backaction noise to manifest itself in- 
creases with increasing signal bandwidth (measurement 
noise hides the backaction), increases with particle mass 
(scenario is less "quantum"), and increases with mea- 
surement inaccuracy (need a good measurement to get 
to the SQL). This picture is illustrated by Figure la (for 
fixed D, B, and m), which shows x rras (t) , S rms , and f* 
together with two stochastic realizations of x (t) . 

In the special case where we choose to integrate the 
measurement results over all of a given observation in- 
terval t, B = 1/t and we recover a single measurement 
scenario. We can then derive a basic figure-of-merit for 
continuous measurements — defining S = \f~D as the mea- 
surement inaccuracy and keeping t as the total observa- 
tion time, we find that we need 



f <3"W- 



(12) 



in order for noise in the measurement record ultimately 
to be dominated by real displacements of the particle due 
to backaction diffusion. Alternatively, we can write down 
something like an inference-disturbance relation: 
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This relation asserts that the product of the integrated 
measurement accuracy and the rms displacement of the 
test particle due to backaction-diffusion must exceed a 
lower bound that increases linearly with time. 

We now turn to consider the estimation of weak forces 
acting on the test particle, based on something like 
Kalman filtering of the measurement record. We begin 
by writing down a modified form of the stationary differ- 
ence equations (still explicit), 

T C — 1 

X r+ i - x r =p r 1 — - (£ r - x' r ) , (14) 

m C 

p r+ i -p r = ar H ^= (£ r - x' r ) . (15) 

(TV C/ 



Here we have simply added a uniform force field a that 
acts on the particle at all times E9]. Our task will be to 
detect a^0. 

Based on the measurement results £ r , we want to up- 
date recursive estimators x r and p r for x r and p r , re- 
spectively. Assuming noiseless readout of £ r , the update 
equations read 



x r +i 
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(16) 
(17) 



Note that we do not include a in the update rules, as it 
is assumed to be an unknown quantity. In order to infer 
a from the recursive estimators x r and j5 r , we must focus 
on the behavior of the quantity rj r = £ r — x' r . As we know 
the distribution for £ r , we can write 



Vr 



(^X ' x^, J I 



-K [o, i] , 



(18) 



where M r [0, 1] is a Gaussian deviate with zero mean 
and unit variance. Hence, rj r will generally be the sum 
of a Gaussian deviate and an uncorrelated process (the 
"error-signal" ) that reflects the cumulative inaccuracy of 
our recursive position estimator. 

Note that if there is no external force acting on the 
particle (a = 0) there is no measurement inaccuracy, 
and the recursive estimator x r can be used perfectly to 
subtract the backaction-induced diffusion from the sig- 
nal £ r (see Figure lb). Indeed, demonstrating the ability 
to do so would seem to be the only good way of ver- 
ifying that an actual experimental broadband position 
measurement reaches the "quantum regime" defined by 
equation @— unlike the scenario of optical measure- 
ments of quadrature amplitudes, there is no way directly 
to compare signal and meter beams to determine the fi- 
delity of the measurement. If there is external force, the 
subtraction will not be perfect and a will reveal itself in 
nonstationary behavior of r\ r . 
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Note as we ll that the variance of the residual noise 
process y ~^-N r [0, 1] can in principle go to zero in the 
limit of good measurements. At first thought this might 
seem to have profound implications for the detection of 
weak classical forces, but in fact the quantitative analysis 
below recovers the usual SQL for weak force detection. 

Before proceeding any farther, let us make the transi- 
tion to the continuous limit (r — > 0, a — > oo, D = const) 
jl3| and collect together a complete set of SDE's for the 
physical quantities and statistical estimators: 



dx = — dt, 
m 



d£ = xdt 



dx = —dt, 

m 



dp = adt + \ j ~^j=jdW 7 




dp = — (d£ — xdt) . 



(19) 

(20) 
(21) 



In these terms, d^/dt represents the broadband measur- 
ment record. Note that we could also define an auxilliary 
variable e = x — x, in terms of which 



de = — (p — p) dt, 
m 



dp=-edt+d—dW. 



(22) 
(23) 



E = 



aD 




Setting Y./N = 1 we obtain 
7 



t 



mD . h ' 
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(27) 



(28) 



For fixed m and t, it seems reasonable that the optimal 
choice of D should be simply related to fit 2 jm. If we 
substitute D = r] 2 ht 2 /m into ( |2S| ) and then (numerically) 
minimize over rj, 

a.nin = min J ^ \r) - rj 2 sin (1/?/) 
r) V 2i 

Hence we are able to recover (within a factor of tt) the 
usual SQL for detecting a weak classical force via re- 
peated position measurements on a free mass 

The author would like to thank J. W. Hartman and 
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from the National Science Foundation under Grants No. 
PHY94-07194 and PHY97-22674. 



Our estimator for the external force a will be propor- 
tional to the time integral of drj = d£ — xdt, so what we 
really need to know is the rms behavior of e (t) . Devel- 
oping its SDE, 



de 


= dx 
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(24) 



This shows that e (t) behaves like a harmonic oscillator, 
with mass m and natural frequency u>o = y , sub- 
jected to a constant external force a. So with the initial 
condtion x — x, p — p, we expect 



, . aD 
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(25) 



Getting back to the time integral of dr/, 

f^ )M ' + ^H iw(t,) 

(26) 

from which we identify the signal E and rms noise N as 
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